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We propose to experimentally test the nonclassicality of 
quantum states through homodyne tomography. For single- 
mode states we check violations of inequalities involving the 
photon-number probability. For two-mode states we test 
the nonclassicality by reconstructing some suitable number- 
operator functions. The test can be performed with available 
quantum efficiency of homodyne detection, by measuring the 
pertaining quantities on the corresponding noisy states. 

1998 PACS number(s): 03.65.-w, 03.65.Bz, 42.50.Dv, 42.50-p 



I. INTRODUCTION 

The concept of nonclassical states of light has drawn 
much attention in quantum optics jl] . The custom- 
ary definition of nonclassicality is given in terms of the 
Glauber-Sudarshan P-function: a nonclassical state does 
not admit a regular positive P-function representation, 
namely, it cannot be written as a statistical mixture of 
coherent states. Such states produce effects that have no 
classical analogue. These kinds of states are of funda- 
mental relevance not only for the demonstration of the 
inadequacy of classical description, but also for applica- 
tions, e.g., in the realms of information transmission and 
interferometric measurements [|||||l0|]. 

In this paper we are interested in testing the nonclas- 
sicality of a quantum state by means of an operational 
criterium, which is based on a set of quantities that 
can be measured experimentally with some given level 
of confidence, even in the presence of loss, noise, and 
less-than-unity quantum efficiency The positivity of the 
P-function itself cannot be adopted as a test, since there 
is no method available to measure it. The P-function is a 
Fourier transform on the complex plane of the generating 
function for the normal-ordered moments; hence, in prin- 
ciple, it could be recovered by measuring all the quadra- 
ture components of the field, and subsequently perform- 
ing an (deconvolved) inverse Radon transform |12] . Cur- 
rently, there is a well-established quantitative method for 
such a universal homodyne measurement, and it is usu- 
ally referred to as quantum homodyne tomography (se e 
Rcf. |l3| for a review). However, as proven in Ref. 
only the generalized Wigner functions of order s < l—r/^ 1 



can be measured, r\ being the quantum efficiency of ho- 
modyne detection. Hence, through this technique, all 
functions from s — 1 to s — cannot be recovered, i.e., 
we cannot obtain the P-function and all its smoothed 
convolutions up to the customary Wigner function. For 
the same reason, the nonclassicality parameter proposed 
by Lee || , namely, the maximum s-parameter that pro- 
vides a positive distribution, cannot be experimentally 
measured. 

Among the many manifestations of nonclassical effects, 
one finds squeezing, antibunching, even-odd oscillations 
in the photon-number probability, and negativity of the 
Wigner function p|-p| pO|Jl~5| -p^7[| . Any of these features 
alone, however, does not represent the univocal criterium 
we are looking for. Neither squeezing nor antibunching 
provides a necessary condition for nonclassicality Q| . The 
negativity of the Wigner function, which is well exhibited 
by the Fock states and the Schrodinger-cat-like states, is 
absent for the squeezed states. As for the oscillations 
in the photon-number probability, some even-odd oscil- 
lations can be simply obtained by using a statistical mix- 
ture of coherent states ||l8| . 

Many authors ||7|,p|jl 1| have adopted the nonpos- 
itivity of the phase-averaged P-function F(I) — 
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P(/ 1 / 2 e l< ^) as the definition for a nonclassical 



state, since F(I) < invalidates Mandel's semiclassical 
formula JjJ of photon counting, i.e., it does not allow a 
classical description in terms of a stochastic intensity. 
Of course, some states can exhibit a "weak" nonclas- 
sicality namely, a positive F(I), but with a non- 
positive P-function (a relevant example being a coher- 
ent state undergoing Kerr- type self-phase modulation). 
However, from the point of view of the detection the- 
ory, such "weak" nonclassical states still admit a classical 
description in terms of having the intensity probability 
F(I) > 0. For this reason, we adopt nonpositivity of 
F(I) as the definition of nonclassicality. 



II. SINGLE-MODE NONCLASSICALITY 

The authors of Refs. have recognized relations 

between F(I) and generalized moments of the photon 
distribution, which, in turn, can be used to test the non- 
classicality. The problem is reduced to an infinite set 
of inequalities that provide both necessary and sufficient 
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conditions for nonclassicality M . In terms of the photon- 
number probability p(n) — (n\p\n) of the state with den- 
sity matrix p, the simplest sufficient condition involves 
the following three-point relation for p(n) 

B(n) = (n + 2)p(n)p{n + 2) 

-(n + l)[p(n+ l)] 2 < . (1) 

Higher-order sufficient conditions involve five-, seven-, 
. . . , (2k + l)-point relations, always for adjacent values 
of n. It is sufficient that just one of these inequalities be 
satisfied in order to assure the negativity of F(I). Notice 
that for a coherent state B(n) = identically for all n. 

In the following we show that quantum tomography 
can be used as a powerful tool for performing the non- 
classicality test in Eq. (Q). For less-than-unity quantum 
efficiency (rj < 1), we rely on the concept of a "noisy 
state" Qr) , wherein the effect of quantum efficiency is as- 
cribed to the quantum state itself rather than to the de- 
tector. In this model, the effect of quantum efficiency 
is treated in a Schrodinger-like picture, with the state 
evolving from g to g„, and with r\ playing the role of a 
time parameter. Such lossy evolution is described by the 
master equation 

p 

d t g(t) = - {2ag(t)a f - a)ag(t) - g(t)a)a) , (2) 

wherein g(t) = g n with t = — lnry/r. 

For the nonclassicality test, reconstruction in terms of 
the noisy state has many advantages over the true-state 
reconstruction. In fact, for nonunit quantum efficiency 
r\ < 1 the tomographic method introduces errors for 
p(n) which are increasingly large versus n, with the addi- 
tional limitation that quantum efficiency must be greater 
than the minimum value r\ = 0.5 On the other 

hand, the reconstruction of the noisy-state probabilities 
Pn(n) — (n\p 71 \n) does not suffer such limitations, and 
even though all quantum features are certainly dimin- 
ished in the noisy-state description, nevertheless the ef- 
fect of nonunity quantum efficiency does not change the 
sign of the P- function, but only rescales it as follows: 

P(z)^P 7l (z) = -P(z/r, 1 / 2 ). (3) 
V 

Hence, the inequality (|l|) still represents a sufficient con- 
dition for nonclassicality when the original probabilities 
p(n) — (n\p\n) are replaced with the noisy-state proba- 
bilities p v (n) — (n\pri\n), the latter being given by the 
Bernoulli convolution 

PnH = jt(^)v n (.l-v) k - n p(k). (4) 

Hence, when referred to the noisy-state probabilities 
p v (n) , the inequality in Eq. (|l|) keeps its form and simply 
rewrites as follows 



B v (n) = (n + 2)p v (n)p v (n + 2) 

-(n+l)[p„(n+l)] 2 <0. (5) 

According to Eq. (|5|), the quantity B n (n) is nonlinear 
in the density matrix. This means that B n (n) cannot be 
measured by averaging a suitable kernel function over the 
homodyne data, as for any other observable Ej|. Hence, 
in the evaluation of B n (n) one needs to tomographically 
reconstruct the photon- number probabilities, using the 
kernel functions 

= 2 J £tf( n _ ) (2, + 1)!^ 

x Re {D_ (2 „ +2) (-2iKx)} , (6) 

where D a (z) denotes the parabolic cylinder function and 
k = ^/ 77/(277 — 1). The true-state probabilities p(n) are 
obtained by averaging the kernel function in Eq. (g) over 
the homodyne data. On the other hand, the noisy-state 
probabilities Prj(n) are obtained by using the kernel func- 
tion in Eq. (|6|) for r\ = 1, namely without recovering the 
convolution effect of nonunit quantum efficiency. Notice 
that the expression (^) does not depend on the phase of 
the quadrature. Hence, the knowledge of the phase of the 
local oscillator in the homodyne detector is not needed 
for the tomographic reconstruction, and it can be left 
fluctuating in a real experiment. 

Regarding the estimation of statistical errors, they are 
generally obtained by dividing the set of homodyne data 
into blocks. However, in the present case, the nonlinear 
dependence on the photon number probability introduces 
a systematic error that is vanishingly small for increas- 
ingly larger sets of data. Therefore, the estimated value 
of B(n) will be obtained from the full set of data, instead 
of averaging the mean value of the different statistical 
blocks. 

In Figs. [l|-|7| we present some numerical results that 
are obtained by a Monte-Carlo simulation of a quantum 
tomography experiment. The nonclassicality criterium 
is tested either on a Schrodinger-cat state \ip(a)) cx 
(|a) + | — a)) or on a squeezed state \a,r) = D(a)S(r)\0), 
wherein \a), D(a), and S(r) denote a coherent state with 
amplitude a, the displacement operator D(a) — e aa _Qa , 

and the squeezing operator S(r) — e 7 '^ ~ a respec- 
tively. Figs. show tomographically-obtained values of 
B(n), with the respective error bars superimposed, along 
with the theoretical values for a Schrodinger-cat state, 
for a phase-squeezed state (r > 0), and for an amplitude- 
squeezed state (r < 0), respectively. For the same set of 
states the results for B n (n) [cf. Eq. (j^)] obtained by to- 
mographic reconstruction of the noisy state are reported 
in Figs. ^-|^. Let us compare the statistical errors that 
affect the two measurements, namely, those of B(n) and 
B^n) on the original and the noisy states, respectively. 
In the first case (Figs. [5] ||) the error increases with n, 
whereas in the second (Figs. ^-^) it remains nearly con- 
stant, albeit with less marked oscillations in B v (n) than 
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FIG. 1. Tomographic measurement of B(n) (dashed 
trace) with the respective error bars (superimposed in 
grey-shade) along with the theoretical values (solid trace) for 
a Schrodinger-cat state with average photon number n = 5. 
The quantum efficiency is r\ = 0.8 and the number of simu- 
lated experimental data used for the reconstruction is 10 7 . 

those in B{n). Fig. ^| shows tomographically-obtained 
values of B n (n) for the phase-squeezed state (cf. Fig. ||), 
but for a lower quantum efficiency rj = 0.4. Notice that, 
in spite of the low quantum efficiency, the nonclassical- 
ity of such a state is still experimentally verifiable, as 
■877(0) < by more than five standard deviations. In 
contrast, for coherent states one obtains small statistical 
fluctuations around zero for all n. We remark that the 
simpler test of checking for antibunching or oscillations in 
the photon-number probability in the case of the phase- 
squeezed state considered here (Figs. ||, ||, and ^) would 
not reveal the nonclassical features of such a state. 




FIG. 2. Tomographic measurement of B(n) (dashed 
trace) with the respective error bars (superimposed in 
grey-shade) along with the theoretical values (solid trace) for 
a phase- squeezed state with n — 5 and n sq = sinh 2 r — 3 
squeezing photons. The quantum efficiency is 77 = 0.8 and 
10 7 simulated experimental data have been used for the re- 
construction. 
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III. TWO-MODE NONCLASSICALITY 

Quantum homodyne tomography can also be employed 
to test the nonclassicality of two-mode states. For a 
two-mode state nonclassicality is defined in terms of non- 
positivity of the following phase-averaged two-mode P- 
function [ffll: 
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#iP(7 1 1/2 e*,/ 2 1/2 e < W 1 +^) . (7) 



In Ref. |TlJ] it is also shown that a sufficient condition for 
nonclassicality is 

C = ((ni - h 2 f) - {{fix - n 2 )f - (ni + n 2 ) < , (8) 

where fix and 77,2 are the photon-number operators of the 
two modes. 

A tomographic test of the inequality in Eq. (||) can 
be performed by averaging the kernel functions for the 
operators in the ensemble averages in Eq. (^) over the 
two-mode homodyne data. For the normal-ordered field 



FIG. 3. Same as in Fig. ^, but for an amplitude-squeezed 
state. 



operators one can use the Richter formula in Ref. 
namely 



H n+m (y/2rjx) 
v /(2r ? )"+™("+ m ) 



(9) 



H n (x) denoting the Hermite polynomial and 4> being the 
phase of the fiels with respect to the local oscillator of the 
homodyne detector. Again, as for the kernel function 
in Eq. (||), the value 77 = 1 is used to reconstruct the 
ensemble averages of the noisy state p v . Notice that for 
n = m Eq. (^) is independent on the phase cf), and hence 
no phase knowledge is needed to reconstruct the ensemble 
averages in Eq. (|J). As an example, we consider the twin- 
beam state at the output of a nondegenerate parametric 
amplifier 
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FIG. 4. Tomographic measurement of B v (n) for a 
Schrodinger-cat state with n — 5, degraded by a quantum 
efficiency r\ — 0.8. The number of simulated experimental 
data is 10 7 . 
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FIG. 5. Tomographic measurement of B v (n) (dashed 
trace) with the respective error bars (superimposed in 
grey-shade) along with the theoretical values (solid trace) for 
a phase-squeezed state, which has n — 5 and n sq = sinh 2 r — 3 
squeezing photons, and which has been degraded by a quan- 
tum efficiency r\ — 0.8. For the reconstruction a sample of 



(10) 10 simulated experimental data have been used. 



where \n) ® \n) denotes the joint eigenvector of the num- 
ber operators of the two modes with equal eigenvalue n, 
and the parameter A is related to the gain G of the ampli- 
fier by the relation |A| 2 = 1 — G _1 . The theoretical value 
of C for the state in Eq. @ is C = -2|A| 2 /(1-|A| 2 ) < 0. 
A tomographic reconstruction of the twin-beam state in 
Eq. ( |l0| ) is particularly facilitated by the self-homodyning 
scheme, as shown in Ref. [p2| . With regard to the effect of 
quantum efficiency r\ < 1, the same argument still holds 
as for the single-mode case: one can evaluate C v for the 
twin-beam state that has been degraded by the effect of 
loss. In this case, the theoretical value of C v is simply 
rescaled to C v = -2r? 2 |A| 2 /(l - |A| 2 ). 

In Fig. |^ we report C n vs. l — r/,r] ranging from 1 to 0.3 
in steps of 0.05, for the twin-beam state in Eq. ( |l0| ) with 
|A| 2 = 0.5, corresponding to the total average photon 
number equal to 2. The values of C n result from a Monte- 
Carlo simulation of a homodyne tomography experiment 
with a sample of 4 x 10 5 data, using the theoretical joint 
homodyne probability of the state \x) 



p v (xi,x 2 , <t>\,<t>i) 



with 



2 exp 



(x 1 +x 2 ) 



(ci 2 +4A 2 )(d 2 _ 2 



\l±z\ 2 



4A 2 ) 



(11) 
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0i and fa denoting the phases of the two modes relative 
to the respective local oscillator. Notice that the nonclas- 
sicality test in terms of the noisy state gives values of C, ; 
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FIG. 6. Same as in Fig. q, but for an amplitude-squeezed 
state. 



that are increasingly near the classically positive region 
for decreasing quantum efficiency rj. However, the statis- 
tical error remains constant and is sufficiently small to 
allow recognition of the nonclassicality of the twin-beam 
state in Eq. (JlOj) up to rj — 0.3. 



IV. CONCLUSIONS 

We have shown that quantum homodyne tomogra- 
phy allows one to perform nonclassicality tests for var- 
ious single- and two- mode radiation states, even when 
the quantum efficiency of homodyne detection is rather 
low. The method involves reconstruction of the photon- 
number probability or of some suitable function of the 
number operators pertaining to the noisy state, namely, 
the state degraded by the less-than-unity quantum ef- 
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FIG. 7. Same as in Fig. || but here for a quantum efficiency 
of r/ — 0.4, and a sample of 5 x 10 7 simulated experimental 
data. 
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FIG. 8. Tomographic measurement of C v as defined in 
Eq. (Eh and modified by the quantum efficiency for the 
twin-beam state in Eq. (101). The respective error bars are 
shown in the grey shade and |A| 2 = 0.5 corresponding to an 
average of 2 total photons. The results are shown for various 
values of the quantum efficiency r\ (in steps of 0.05) and for 
each value of n the number of simulated data is 4 x 10 5 . 

ficiency. The noisy-state reconstruction is affected by 
the statistical errors; however, they are sufficiently small 
that the nonclassicality of the state can be tested even 
for low values of r\. For the cases considered in this pa- 
per, we have shown that the nonclassicality of the states 
can be proven (deviation from classicality by many error 
bars) with 10 5 -10 7 homodyne data. Moreover, since the 
knowledge of the phase of the local oscillator in the homo- 
dyne detector is not needed for the tomographic recon- 
struction, it can be left fluctuating in a real experiment. 
Hence, we conclude that the proposed nonclassicality test 
should be easy to perform experimentally. 
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